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The structure and dynamics of steady state sedimentation of semi-concentrated (<j!> = 0.10) 
monodisperse spheres are studied in liquid fiuidized beds. Laser turbidity and particle imaging 
methods are used to measure the particle velocity fluctuations and the steady state concentration 
profiles. Using a wide range of particle and system sizes, we find that the measured gradients 
V(f), the fluctuation magnitudes (t„, and their spatial correlation lengths ^, are not uniform in the 
columns - they all show strongly 2:— dependent proflles. These profiles also display a scaling in which 
results from different particle sizes collapse together when plotted in the forms —aV(f>[z), £^[z)/a, 
and (T„(z)/up, demonstrating the universality of the particle dynamics and structure in steady state 
sedimentation. Our results are also used to test a recently proposed model for the correlation lengths 
(^{z) in terms of the concentration stratification V4>{z) [P.J. Mucha and M.P. Brenner, Phys. Fluids 
15, 1305 (2003)], e,{z) = coa[ii)S{(l>)f'^[-aV(l>{z)\-'^''\ We find that the correlation lengths pre- 
dicted by this model are in very good agreement with our measured values, showing that the origin 
of the fluctuation length ^ lies with the concentration stratiflcation V^. 

PACS numbers: 47.15. Gf, 05.40.-a, 47.55.Kf, 82.70.Dd 



I. INTRODUCTION 

The slow sedimentation of a collection of non- 
Brownian spheres in liquids is a fundamental problem 
in physics and is of importance in numerous chemical 
Industrie^. Despite its apparent simplicity, the com- 
plexities of hydrodynamic i nteracti ons between particles 
are still vigorously debatecP'^'^'^. Much of the debate 
surrounds the origin of the characteristic size of veloc- 
ity fluctuations, ^ ~ llacj)^^/^, found in sedimentation 
experiment^^Il] (jjere a is the particle radius and 4> is the 
volume fraction). This lengthscale is typically smaller 
than the cell dimensions and plays a central role in sed- 
imentation dynamics. The magnitude of the particle 
velocity fluctuations (T„ depends directly on ^. In di- 
lute suspensions, Poisson statistics dictates that in a re- 
gion of size ^, the typical concentration fluctuation is 
C0 — \J ^o? its velocity is determined by balanc- 
ing its buoyant weight, a^^pg^ , with its Stokes drag, 
Girrj^ay, giving ay cx y/cjj^/a 'J ^ Similarly, particle diffu- 
sion coefficients, D ~ S^ay, also depend directly on 
Due to its central role, an understanding of the origin of 
the ffuctuation correlation length is essential for even a 
basic description of sedimentation. 

Recent theoriej ^° * ^^ * ^^ * ^'^ ^ have taken a new ap- 
proach to understanding the lengthscale of velocity ffuc- 
tuations by considering the possible effects of a small 
stratificatio n in m ean column concentration, Vcj). Mucha 
and Brennei'i^'^ argue that because velocity ffuctuations 
of size I are driven by buoyancy forces acting on local 
concentration ffuctuations cr^, when there is a concen- 
tration variation across the ffuctuating region due to the 
stratific ation that is of order the fiuctuation itself, 
= y^c/ja^/P, then the buoyancy mismatch of the fluc- 
tuating region is reduced, and motions are suppressed. 
The largest region for which motion occurs then is of 



order I — —a^/Vf/), or 

l/a^co[<pSm'^'[-aV<j>r'^', (1) 

which relates the correlation size / to the concentra- 
tion stratiflcation V(/). Here, the static structure factor 
S^ip) accounts for excluded volume affects in non-dilute 
sampled^. Results in support of this modeP^, using nu- 
merical simulations, have focussed on the diffusion coef- 
flcients at the transition region near the top particle in- 
terface. The generality of the Mucha and Brenner (MB) 
model, however, suggests it is not conflned solely to inter- 
facial properties, and may provide a complete description 
for the fundamental length ^ over the entire sedimenting 
column. 

The most direct test of the MB model, as performed 
in this article, would be a direct comparison between 
measured fluctuation lengths and the calculated strat- 
iflcation lengths /, over all heights in a sedimenting par- 
ticle column. This test has yet to be done, principally 
due to a lack of data for both ^(z) and V(/)(z). More- 
over, the few available data for V0 far from the interface, 
usin g si mulations, yield seemingly contradictory results; 
LadcpS found negligibly small gradients, but Mucha and 
Brennei^^l found much larger ones. We note that since 
Mucha and Brenner considered larger systems than Ladd, 
this contradiction may be resolved if gradients have an 
as yet unidentifled system size dependency. 

In this article we examine steady state sedimenta- 
tion using a wide range of particle and system sizes. The 
samples have the same column height, H, and the same 
mean concentration, c/jq = 10.0%, but different particle 
sizes a, which range over a factor of ~ 4.5 so that the 
system size N (cx a^) varies by almost two orders of mag- 
nitude. We flrst show that all of the samples show some 
degree of stratiflcation. The magnitude of the stratiflca- 
tion increases with (normalized) cell size until its aver- 
age becomes of order V(/> ~ (po/H. We then show that 




FIG. 1: Photograph of a random sampHng of a — 154.5 fim 
spheres resting on a glass slide. 



there is a scaling regime in which the (height-dependent) 
concentration gradients, the velocity correlation lengths, 
and the velocity fluctuation magnitudes all collapse onto 
master curves when scaled as —aV4'{z), ^{z)/a, and 
(jy{z)/vp. All of our data are consistent with the MB 
model. We find that the stratification lengths l{z) calcu- 
lated from Eq. ([l]) agree with our measured values ^(z) 
using a single fit constant, cq = 0.82, for three differ- 
ent samples in the scaling regime, and a slightly lower 
value, Co = 0.70, in the non-scaling regime. These find- 
ings validate the physical picture behind the MB model, 
and show that the origin of the fluctuation length ^ lies 
with the concentration stratification Vt/f. 

Our choice of an experimental arrangement to study 
Vcj) is guided by the recent observatio ns of pe rsistently 
evolving gradients during sediment atiorfi2lISE3 xhe evo- 
lution occurs because, starting from a well mixed sam- 
ple (V0 = 0), it takes some time for the stratification 
(V(/) ^ 0) to become established. During this evolu- 
tion, the column is falling, and often steady state is never 
reached before all the particles lie on the bottonJi^l, Xo 
ensure we reach a steady state, we conduct experiments 
in a fiuidized bed apparatus, where liquid is pumped up- 
wards to balance the settling of the particleJ^^E^l xhe 
fiow balancing ensures that the particle column never col- 
lapses, and a steady state is eventually established. 



II. EXPERIMENT DESCRIPTION 

A. Particles and Fluids 

We use spherical glass beads of mean radii a — 
54.5, 103.5, 154.5, and 245 /im. The beads have been spe- 
cially filtered by the manufacturer (Mo-Sci corp.) to be 
of extremely low polydispersity in size. To demonstrate 
this, we show in Fig. [l]a micrograph of a collection of the 
a — 154.5 beads. The particles lay on a flat plate and 
through a gentle shaking quickly organize into highly or- 
dered forms, indicative of low levels of polydispersity aa- 
To accurately determine CTq, for each set of particles we 




10 20 30 40 50 
Particle Number / 



FIG. 2: The settling velocities of 50 isolated spheres falling in 
a water/glycerol mixture. The solid line is the average settling 
rate, v — 2.218 mm/s, and the dashed lines the standard 
deviation, = 0.047 mm/s. The polydispersity in velocity is 
(Jv/v = 2.13%, corresponding to a particle size polydispersity 
of aa/a = 1.07%. 



TABLE I: Particle, fluid, and fiuidized bed properties, a is the 
average particle radius, aa/a the measured size polydispersity, 
■q the pure solvent viscosity, vo the Stokes settling velocity, 
Re the particle Reynolds number, H the total height of the 
fiuidized particle column, and (po the average volume fraction 
of the fiuidized bed. 

a{ij.m) aa/a rj (cp) vo{mm/s) Re H (cm) (po 



54.5 0.016 17 0.548 

103.5 0.015 27 1.21 

154.5 0.012 54 1.34 

245 0.011 87 2.10 



0.002 18.5 0.10 

0.005 18.5 0.10 

0.004 18.5 0.10 

0.006 18.5 0.10 



measure the variation in settling velocities of 50 individ- 
ual spheres. Fig. |2] shows the results for the a = 245 
spheres. The standard deviation of the individual parti- 
cle settling rates is indeed extremely low, ay/v — 2.13%. 
Using the Stokes formulcP', vg = {2/9) A pga'^ /rj, and as- 
suming the particles are all of equal density, the variation 
in settling rates correspond to a variation in particle size 
of aa/a = 1.07 %. This procedure is repeated for all 

4 particle sizes, with the results listed in Table I. In all 
cases, the measured value of the size polydispersity aa/a 
is less than 1.7%. 

The glass beads are dispersed in viscous solutions of 
glycerol and water chosen so that inertial forces during 
fluidization are negligible. The particle Reynolds num- 
bers, listed in Table I, are all of order Re = voap/r] « 

5 X 10"'' <^ 1. Additionally, particle motions always oc- 
cur at very high Peclet numbers (Pe ~ 10^), so that 
Brownian diffusion is negligible. The temperature is at 
the ambient value T = 21°C. 




FIG. 3: Photographs of a fluidized particle bed comprised of 
D = 490 /xm spheres. LHS - whole cell. RHS - closeup of the 
top of the column. 



B. Fluidized Bed 

The fluidized bed, shown in Fig. |3] consists of a fluid 
and particle flUed glass cell at the bottom of which fluid is 
pumped upwards to counteract the particle settling and 
fluidize the particles. The sample cell is a rectangular 
glass tube of dimensions TxVl^xi7=:8x 80 x 305 
mm. The overflow liquid at the top of the cell recirculates 
back into the pump, forming a closed loop. To enable a 
uniform flow into the cell, a 2 cm thick nylon mesh is 
packed with 0.5 mm diameter beads and glued across 
the entrance to the cell at the bottom. 

In the experiments described below, we use this setup 
to fluidize four different particle size samples, as listed in 
Table I. To examine the speciflc effects of particle size, 
all four samples share the same average particle concen- 
tration, (j)o = 0.10, and same total height of the particle 
column H. With the pump off, the spheres in all 4 sam- 
ples form a sediment at the bottom ~ 2.9 cm tall. When 
the pump is on, the particles expand upward, flUing a re- 
gion above the bottom up to a height dependent upon the 
pumped fluid velocity Vp (see Fig. For all 4 samples, 
we set Vp to a value that expands the particle column to 
a total height H ~ 18.5 cm, so that the average volume 
fraction is approximately the same for all of the samples, 
00 = (IJB.CP * (2.9/18.5) - 0.10, where ^^cp = 0.638 is 
the assumed value for the sediment at the random close 
packing density. 



C. PIV Imaging System: Velocity Flow Maps 

Particle velocities are measured using the technique 
of particle image velocimetry (PTVj^. The apparatus 
consists of a (1008 x 1024 pixels) CCD camera, a syn- 
chronized stroboscope illuminating the cell from behind, 
and specialized image processing hardware and software 
from Dantec Instruments. The depth of fleld of the cam- 
era lens is ~ 5 mm. A large cross section of the cell is 
imaged so that several thousand particles can be simul- 
taneously studied. Velocity maps consisting of 62 x 62 
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FIG. 4: Velocity vector maps of a stable fluidized bed of av- 
erage concentration (po — 0.100. Figures 1 (a) and (b) corre- 
spond to respective positions near to the top and near to the 
middle of the particle column. The single arrow on the right 
gives the corresponding scale of the Stokes settling velocity 
vq. Note that the velocity scale in (a) is magnified relative to 
(b) by a factor of 2 for clarity. 



vectors are extracted by comparing two closely timed pic- 
tures using standard PIV techniques. Each vector is the 
average velocity of two to four spheres. 

Figure [4] shows typical velocity vector maps from a 
stable fluidized bed, where (a) corresponds to a position 
near the top and (b) to a position near the middle of the 
particle column. For scale, we also show the magnitude 
of the Stokes settling velocity vq. Both velocity maps 
show regions moving upwards and downwards, and, as in 
all of our samples, the magnitudes of the velocities are 
significantly larger near the middle than near the top. 

To quantify these observations, we measure the ve- 
locity maps at different heights z along the particle col- 
umn. We extract from the velocity maps the mean ve- 
locities, Vz = {vi,z), and the root mean square (rms) 
velocity fluctuations, tr^ = ((wi,z — Vz)^)^^^, where (. . .) 
represents an ensemble average over ~ 50 vector maps of 
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FIG. 5: Correlation functions of the velocity fluctuations C{z) 
as a function of distance z, near mid-height in the a — 54.5 /im 
particle fluidized bed. The inset shows the logarithm of C'{z), 
the solid line being a fit to C{z) = exp(— (jz/^)^'^), from which 
we extract values of the correlation lengths ^. 



3844 vectors each. 

The typical velocity vector maps in Fig. |4] also dis- 
play large regions where the velocity vectors are spatially 
correlated. To quantify this, we calculate the normalized 
spatial correlation functions of the vertical velocity Vz, 
defined as C{z) = {v z{Q)v z{z)) / {v ^{QY) ■ Figure [s] shows 
typical results for C(z), measured near mid-height in a 
fluidized bed of a = 54.5 ^m particles. The inset shows 
that the empirical form C{z) = exp (— (z/^)^-^) fits the 
data well, enabling us to extract values for ^, the charac- 
teristic (vertical) correlation length of the velocity fluc- 
tuations. 



1. Fluidized Bed Stabihty 

Before using the fluidized bed apparatus for data col- 
lection we need to ensure that the bed is stable in time 
and that there are no large scale convective currents in- 
dicative of an unbalanced (i.e. tilted) cell. To do this we 
use the PIV imaging method and examine particle ve- 
locities over different regions in the cell. Figure |6] shows 
typical results for v and cr„, both perpendicular and par- 
allel to gravity, from PIV vector maps spanning a 4 x 4 
cm imaging window near mid-height. All quantities ran- 
domly fluctuate, but they also maintain time-invariant 
mean values indicative of a stable system. Additionally, 
both velocity averages fluctuate about v — 0, explicitly 
demonstrating that particle sedimentation downwards is 
balanced by the fluidizing counter flow upwards. During 
the experiments, we monitor the mean velocities to en- 
sure stability, and we also record the position of the top 
particle/fluid interface, clearly visible in Fig. [s] to ensure 
that it does not rise or fall in the cell over time. 

As a way of further demonstrating the stability of 
the bed, we apply an uncontrolled shock to the system 




-0.2 - 



t/r. 



FIG. 6: Fluidized bed stability: Average velocities v, and 
velocity fluctuations (t„, parallel (thick lines) and perpen- 
dicular (thin lines) to gravity as a function of time t/rst, 
where Tst = a/vo is the Stokes time. The imaging window 
is near mid-height in a fluidized bed of particles of mean size 
a — 154.5 nm. 
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FIG. 7: Response of a stable fluidized state to a large per- 
turbation. The velocity fluctuations cr„, both parallel and 
perpendicular to gravity, as a function of time t/rst- At time 
t = a rapid and large scale impulse is applied throughout 
the particle bed to disrupt the particles. 



and examine whether or not it returns back to its original 
state. Figure [T] shows the test results in which we inject 
a large {'^ 3 cm radius) air bubble into the middle of the 
particle bed to cause a rapid and large scale disruption 
to the particles. The particle velocity fluctuations a^, 
from a 4 X 4 cm subsection of the cell at mid-height, are 
driven to values ai t — that are more than 10 times 
their values before the shock. Nevertheless, this distur- 
bance rapidly decays away until, for times t > 250x^4 , the 
system returns to its prior stable state. 



D. Light Scattering System; Particle 
Concentrations. 



Local particle volume fractions are determined from 
the local optical turbidity through the particle column. 
We pass an expanded He-Ne laser beam, of diameter 
~ 0.5 cm, through the fluidized bed at a particular height 
z, and measure the transmitted laser intensity It{z) us- 
ing a CCD camera. Results for a typical sample are 
shown in Fig. [sj Figure [sj^a) shows the measured in- 
tensity profiles on the CCD at different heights in the 
particle column. It's evident that the transmitted inten- 
sity varies with height, with the highest intensity at the 
top, and the lowest at the bottom. To find the corre- 
sponding particle concentrations, we need a calibration 
reference for the dependence of the transmitted intensity 
on particle concentration. To do this, we make several 
reference fluidized beds of differing average concentra- 
tion (j) and measure the transmitted intensity patterns 
at mid-height. The concentration values are calculated 
from the ratio of the heights of the sedimented particles, 
at Hsed and ((> = 0.638, to the height of the fluidized col- 
umn, H , as 4> = 0.638* {Hsed/H). The results are shown 
in Fig. |8jb). No efforts are made to model the observed 
dependence of / on 0, our aim being simply to have a cal- 
ibration reference. Fits of I{z) and I{4>) in Fig. [sj^a) and 
(b) to Gaussian functions yield peak intensity values Ipk , 
plotted in Fig. |8jc). By comparing the measured values 
of Ipk{z) in our fluidized bed with the reference values 
of Ipk (0) we are able to determine the height dependent 
concentrations in the fluidized beds, as shown in 

Fig. ^d). 



III. RESULTS 



In this section we describe our main results for the 
particle dynamics of 4 fluidized beds that all share the 
same average particle concentration, (f>Q = 0.10, and col- 
umn height, H = 18.5 cm. The beds differ, however, 
in the size of the fluidized particles, which range from 
54.5 < a < 245 /jva. As a consequence, the beds differ in 
two significant ways. First, the normalized bed heights, 
H* = H/a, range from 755 for the largest particles to 
3364 for the smallest. Second, the total number of par- 
ticles fluidized, N = [{3/AT:)<j)oTWH] x a^^ range from 
192, 000 < < 17, 500, 000, as listed in Table II. A com- 
parison of the results for different particle sizes will allow 
us to examine the specific dependencies of fiuidization on 
the system size N, and the effective column height H* . 
For each fluidized bed we present results for the height de- 
pendence of the particle concentration ^(z), the velocity 
fluctuations cr^(z), and the velocity correlation lengths 
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FIG. 8: (a) Transmitted laser intensity I vs. height z through 
a fluidized bed of radii a — 103.5 nm particles at average 
volume fraction 00 = 0.10. The solid lines are fits to Gaussian 
functions, (b) Intensity vs. (p calibration. Transmitted laser 
intensities measured at mid-height in fluidized beds of average 
volume fraction 0.046 < (0) < 0.142 as labelled. For clarity, 
only every 25'th point is plotted. The solid lines are fits to 
Gaussian functions, (c) Peak transmitted intensity Ipk{z) vs. 
height z as extracted from Gaussian fits to the data in (a). 
Inset: Calibration curve of peak transmitted intensity Ipk{<t>) 
vs. cj> as extracted from Gaussian fits to the data in part 
(b). (d) Volume fraction vs. height z as determined by 
comparing the I[z) curve in (c) with the calibration curve 
in the inset of (c). 



TABLE II: Measured properties of the fluidized beds de- 
scribed in Table I. a is the particle radius and A*' the total 
number of particles fluidized. Vp/vo is the normalized fluidiz- 



ing pump velocity upwards. 



"/uo is the initial sedimen- 



tation velocity of the top interface, measured just after the 
fluid pump is turned off. 0'°'' is the particle concentration 
measured just below the top interface. The particle concen- 
trations <f)^°^ and 0*°^ are calculated from Vp/vo and i]^°Z/vo 
using the Richardson-Zaki equation v/vo = (1 — 0) = '^. 
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245 
154.5 
103.5 

54.5 


755 
1194 
1769 
3364 


1.92 ■ 10= 
7.67- 10^ 
2.55- 10^ 
1.75- 10^ 


0.605 0.600 
0.659 0.654 
0.727 0.686 
0.797 0.776 


8.4% 
7.0 
6.0 
3.7 


8.7 % 8.9 % 
7.3 7.4 
5.6 6.6 

4.0 % 4.5 



A. Concentration Profiles (^(z). 

We begin by examining the concentration profiles 
0(z) of the four different particle size fluidized beds. The 
laser scattering method, described in section II. D above, 
allows for the accurate determination of the local particle 
volume fraction 0(z) as a function of height z. To do this, 
we flrst fluidize each of the four different particle sized 
samples to the same mean concentration, 0o — 0.100, 
and the same total height, H — 18.5 cm. After waiting 
for transients to die out (a time > H/vq ~ 5 min) we 
obtain the time invariant concentration profiles shown in 
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FIG. 10: The concentration just below the top interface 0'°'', 
o, vs. normalized column height H* , from Fig. [o] Also shown 
are the volume fractions A calculated from measurements of 



the settling velocities of the top interface, 



,,'OP 



and the 



fluidization pump velocity Vp/vo, A, using the RZ equation 
v{(j})/vo = —(1 — (j))^'^ ■ Inset: The initial settling velocity of 



the top interface, 



top 



, and the fluidization pump velocity 



Vp/vo, A, as a function of column height H* . 



Perhaps the most striking result in Fig. [9] is that 
the concentration profiles and the degree of stratification 
change with normalized bed height H* = H/a. That is, 
despite the fact that all of the fluidized beds share the 
same average concentration, (pQ — 0.100, and the same 
total height H = 18.5 cm, there is no universal curve for 
the concentration profiles. Rather, the profiles are found 
to markedly change with normalized height H*. The 
main trend in evidence, as described further below, is for 
the degree of stratification to increase with increasing cell 
height H* , or system size N. 



Interfacial Concentrations. 



FIG. 9: Steady state concentration profiles 1^(2). The 
particle radii, normalized column height, and number of 
particles are (a) a = 245 /xm, H* = 755, A*' = 
1.92B5, (b) 154.5, 1194, 7.67_E;5, (c) 103.5, 1769, 2.55J56, and 
(d) 54.5, 3364, 1.75i?7. The dashed lines show the mean con- 
centrations, (j)o — 0.100. Solid lines are fits to the empir- 
ical form (f>{z) — bo — biz — 626*^^, where 6o,bi,fe2 and 63 
are adjustable fit constants. Insets: concentration gradients, 
-V0 = 61 (6263)e''3^ 



Figs. |9];a-d). 

It is immediately apparent that none of the fluidized 
beds are perfectly homogeneous in concentration. Each 
concentration profile shows the same trend for the con- 
centration to decrease with height, i.e. they are all strat- 
ified. While most of the concentration variation occurs 
in the highest portions of the columns, measurable and 
significant gradients occur at all heights. 



Before examining the details of the gradients in con- 
centration V(j){z) from Fig. |9] we first find it instruc- 
tive to consider a simpler measure of the degree of non- 
uniformity in the beds by considering a measure of the 
variation of in each column. The smallest values of 
are always recorded at the tops of each column, just 
below the sharp interfaces (see Fig. [3] and [Tl]) . To within 
our resolution of 0.5 cm, there is a near discontinuity 
as the concentration drops from a finite value to zero over 
0.5 cm. (Note that much more diffuse int erfaces occur 
for more dilute and/or polydisperse systemd^^E^Ml) ^ •\Ye 
therefore define 0*°^' as the value of (j>{z) measured closest 
to, and just below, the top interface. 

Figure 10 shows our measured values of 0*°'' vs. the 
normalized bed height H* . The overall trend is clear, 
the taller the column H* , the smaller the concentration 
and the greater the degree of stratification. For our 
3 shortest columns, there is a near linear decrease of (j>*°P 
with H* . This linear trend, of course, cannot continue 
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FIG. 11: Photographs of the top portion of a particle column, 
comprised of a = 154.5 fim spheres, taken at different time 
intervals t after having turned off the fluidizing pump. 



indefinitely, as (f> ^ cannot go below zero. Indeed, in 
our tallest bed, the result deviates significantly from the 
initial linear trend. 

To help corroborate the concentration values 0*°^ 
just below the top interface, we also conduct a second se- 
ries of measurements. Because the hindered settling func- 
tion of a collection of spheres is strongly (/)-dependent (the 
Richardson-Zaki (RZ) equatioiP2!is v{4>)/vq — (1 — 0)^-^), 
by measuring the settling velocity of the top interface, 
w*°^/wo, we can independently extract out values for (/>*°''. 
To test this, we first establish a stable fluidized bed, then 
turn off the fluidizing pump and measure the initi al s et- 
tling rate of the top interface, — Wsed/''^o- Figure [ll] il- 
lustrates how these experiments are done. Results for 
— t;*°2/wo vs. H* are shown in the inset of Fig. [lo] The 
corresponding (RZ calculated) values of 0*°^ are shown 
in the main portion of the Figure, and agree to within 
the error bars with the turbidity measurements. In this 
way, we have a second, independent, way of determining 

that further corroborates the non-universality of the 
concentration profiles with particle size a. 



B. Concentration Gradients V(f)[z). 

The insets in Fig. [9](a-d) show the concentration gra- 
dients — V0, as derived from fits to (f){z). The gradients 
are all nearly constant up to mid- height, and increase 
markedly thereafter. They also increase in magnitude 
with decreasing particle size a. To examine this in closer 
detail, we plot in Fig. 12 the scaled values — aV0 vs. 
height z. Remarkably, data for the 3 largest particle 
sizes collapse closely together, showing that in this regime 
there is a new scaling relation, — V0 oc 1/a. This scaling 
also shows that the stratification increases with system 
size. Since the total number of particles in our samples, 
N cx we can equivalently recast the scaling form 

as —W4> oc 7V^/^. Significantly, this form predicts that 
in limit of small system sizes, N Q, the stratifica- 
tion vanishes so that the sample becomes homogeneous 
in concentration. 

It is also readily apparent that the scaling form 
— V(/) cx 1/a predicts an un-physical divergence of — V0 
in the small particle limit a — > 0. Consequently, the scal- 
ing form — aV(/) must eventually break down for small 
enough particles, and indeed, we do see large a large de- 
viation from scaling in Fig. [12] for our smallest particle 
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FIG. 12: Scaling of the concentration gradients, —aVcj) vs. 
height z, from fits to the <f>{z) profiles shown in Fig. [9] 
(— aV</> = ahi -I- {ab2b^)e^^''). The smallest particle size de- 
viates from scahng because the gradient cannot exceed its 
maximum value (see text). The different symbol shapes cor- 
respond to differing particle sizes as defined in Fig. [9] 



size, a — 54.5 /im. To help explain this, we note that 
the largest (linearly averaged over all heights z) value of 
— V(/) that is physically possible, for a bed of height H 
and average concentration ^q, is {—^4')max '^4>o/H. 
Our height averaged gradients in the scaling regime, 
- (V<?!)(z)) = jj j"\/(j3{z)dz - 4.6 • 10~Va, will reach 
this maximum gradient when the particles become small 
enough, of order a* ^ 4.6 • lO~^7?/20o ^ 43 /im. By 
these arguments, the —aW(j) data collapse cannot be 
maintained for bead sizes a < a* = 43 /im, because 
this would require (— V0) to exceed its maximum value 
2(/)o/-ff. Consequently, for beads close to or smaller than 
a* , plots of the product —a\/(j> should deviate below those 
seen in the scaling regime, exactly as we see in Fig. [12] 
for our only set of beads that are close in size to a* . 



1. Stratification Prediction Resulting from Particle Size 
Polydispersity 

Before continuing further, we critically examine the 
degree to which our results for \/(j)(z) may be influenced 
by the non-zero particle size polydispersity aa/a of our 
beads. In general, a concentration stratiflcation can arise 
due to a combination of (1) an inherent physical process 
found in monodisperse beads, and (2) purely the effects 
of particle size polydispersity. Several relevant points 
can be made for our systems. Firstly, we note that all 
4 particle sizes used have similar levels of polydispersity, 
with 1.1 < aa/a < 1.6%, yet the degree of stratifica- 
tion varies greatly between the samples, suggesting that 
polydispersity is not the root cause. Additionally, recent 
computer simulations^'^ compared the sedimentation dy- 
namics and structure of particles with polydispersities 
ranging from 0% (perfectly monodisperse) up to 10%. 
While differences between the monodisperse and 10% 
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FIG. 13: Normalized velocity fluctuations a^, large symbols, 
and the mean particle velocity v, small symbols, parallel to 
gravity as a function of height z in stable fluidized beds. The 
different symbols correspond to the different particle sizes as 
defined in Fig. [9] The solid lines are guides to the eye. 




cases were found, a polydispersity of 2%, slightly larger 
than our samples, yielded results that were, to within 
errors, identical to the monodisperse case. These argu- 
ments suggest that particle size polydispersity is not the 
dominant cause of the observed concentration gradients. 



C. Velocity Fluctuations (t„. 

To measure the particle velocity fields, we use PIV 
methods. Figure [13] displays the values of the particle 
velocities for the 4 different particle size fluidized beds. 
Measurements are taken at multiple heights from near to 
the bottom up to the top, at a height oi H ^ 18.5 cm. 
In all 4 samples the average velocities ?J ~ 0, indicative 
of stable fluidization. Results are also shown for the nor- 
malized velocity fluctuations, Uy/vp, for the four different 
particle size beds. 

In all samples, the fluctuation amplitudes decrease 
with height z in the column. At the top, fluctuations 
become vanishingly small, i.e. the particles are barely 
moving, which we also confirm by direct observation. 

The behavior of the velocity fluctuations also mirrors 
that of the concentration gradients. That is, as found for 

shows (i) a scaling regime, for 



-aV</> in Fig. [12) Fig. [13 

a 3> a*, in which the fluctuation values ay(z)/vp collapse 
together, and (ii) deviations from scaling for a ~ a*. 



D. Velocity Fields-Correlation Lengths 

We next describe the characteristic lengthscale ^ 
of the velocity flucutuations. From the velocity maps 
as measured by PIV, we calculate the spatial corre- 
lation function^ of the vertical velocity fluctuations, 
C{z) = {v^{0)v^{z))/{v^{Of), fit them to the form 



FIG. 14: Velocity correlation lengths plotted as ^ vs. height z 
in (a), and ^/a in (b). The solid lines in (a) are the predictions 
of MB model for I, Eq. (1), with fit constants, in order of 
decreasing particle size, co — 0.82,0.82,0.82 and 0.70. (b) 
The MB predictions for I /a. 



C{z) — exp [(— z/^)^-^], and extract values for for ^ for 
the 4 different particle size fluidized beds are shown in 
Fig.~ 



14 ». 



The correlation lengths are not uniform in height, 
and exhibit a slight decrease towards the top part of the 
column^^. Moreover, at every height z, ^ increases in 
magnitude with increasing particle size. To examine this 
further, we plot in Fig. [l4[^b) the values of the correlation 
lengths normalized by their particle radii, i.e. ^/a. Sig- 
nificantly, the correlation lengths show the same behav- 
ior as the concentration gradients and fluctuation magni- 
tudes. That is, (i) a scaling regime, for a 3> a*, in which 
the correlation lengths (,{z)/a collapse together, and (ii) 
deviations from scaling for a ^ a* . 



E. Connections between Fluidization and 
Sedimentation 

While in this article we focus entirely on the dynam- 
ical properties of fluidized particle beds, it is also of in- 
terest to understand the similarities between fluidization 
and sedimentation. In the lab frame, fluidization and 
sedimentation differ in that fluid is pumped upwards in 
the former, and not pumped at all in the latter. In the 
particle frame of reference, however, the particles in both 
cases are moving through a background fluid at a rate de- 
termined by their buoyant mass and concentration. It is 
not evident, a priori, how the particle dynamics differ 
between fluidization and sedimentation. 
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FIG. 15: Comparison of fluidization and sedimentation, (a) 
the velocity fluctuations a„ and mean velocity w as a function 
of time t. (b) The upward fluid velocity Vp of the fluidized 
bed pump. For times t < 0, the pump is on, and the particles 
are stably fluidized. At time t = 0, the fluidization pump 
is turned off, and the particle sediment downwards until, at 
t = 10 s, the pump is turned back on to re-fluidize the system. 
The fluctuation magnitudes (t„ show no discernible difference 
between fluidization and sedimentation. 



For a direct comparison on a typical sample, we show 
in Fig. [15] results for particle velocity fluctuations for the 
same sample in both fluidized and sedimenting states. 
Fluidization occurs for times leading up to t = sec, 
when the fluid pump is on, and the average particle ve- 
locities in the vertical direction v are very close to zero. 
Sedimentation begins at i = sec, when the fluid pump 
is temporarily turned off, and the particles move down- 
wards with an average velocity ~ —0.9 mm/s. (The 
rapid response of the particles to the turning on and off 
of the fluid pump reflects the lack of particle inertia.) Fi- 
nally, at i = 10 s, after the particles have fallen ~ 1 cm 
(much less than the initial bed height _ff ^ 18 cm), the 
fluidizing pump is turned back on, and the particles are 
re-fluidized. 



The remarkable feature seen in Fig. 15 is that the 
fluctuation magnitudes show no discernible difference 
throughout the transition from fluidization to sedimen- 
tation and back to fluidization. We interpret this as ev- 
idence that the particle fluctuation dynamics during flu- 
idization and sedimentation are very similar. To examine 
this further, in the following two subsections we compare 
our fluidization results for ^ and Uy to published results 
for the same quantities measured in sedimentation exper- 
iments. 



1. Comparison to ^ from Sedimentation Experiments 

Results for the correlation lengths £,/a in sediment- 
ing samples, covering a wide range of volume fractions 
0.0001 < (j) < 0.50, have been published by Segre et 
dlSM^ The methodology by which ^ was measured differs 
from that used here during fluidization. In the sedimen- 
tation experiments, because of the limited time during 
which the particle column is falling, it was not possible 
to measure the fluctuations at different heights. Rather, 
a fixed camera recorded data continuously while the col- 
umn fell past. This produced a single value of ^ that 
represents correlation values averaged over a large por- 
tion of the column, 0.25 < z/H < 0.75. When results 
from samples of various average concentrations 00 were 
collected, the data were consistent with the scaling rela- 
tion ^5 = lla^o ^^^J^ For a meaningful comparison to 
our fluidization results, we thus need to first perform a 
sim ilar height-averaging procedure to our data from Fig. 
[l4| i.e. (^(2)) — jj Jq 25h ^(-^)^-^' "^it^ results plotted 



m the form {£,{z)) vs 
data fit well to {^{z)) 



16 The 



a in the inset of Fig. 
(25.8 ± 2)a, in remarkably good 
agreement with the scaling relation from sedimentation. 
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23.7a. 



We emphasize again that the scaling relation found 
in sedimentation, ^5 — lla(j)Q represented a single, 
height-averaged, correlation length for a column of mean 
concentration (pQ. By comparison, our fiuidization results 
in Figs. |9]and 14 display the complete height dependence 
of both ^ and (f) from within a single column. By com- 
bining these data, we can ask a new question. Can the 
scahng form ^ = lla(j)~^^^ be used to predict the local 
values of ^ from the local values of </> within a single col- 
umn? To answer this, we plot in Fig. [I6] results from 
all four fiuidized beds for ^{z)/a vs. the corresponding 
local volume fractions <j}{z). It is immediately apparent 
that the scaling form does not describe the data, and has 
the opposite dependence on (p. Moreover the data points 
themselves do not collapse together, suggesting that, un- 
like the column averaged values {£,{z)), the values of ^ 
within a column are not a unique function of (p alone. 



2. Comparison to from Sedimentation Experiments. 

Results for the velocity fiuctuations Uy in semi-dilute 
sedimenting samples have been published by numerous 
groups. Guazzelli24- was the first, and found for a sample 
at the same concentration used here, pQ — 0.10, that 
f^vhsed = 0.97±0.19. Later, Segr^ measured (Ty/vsed = 
0.70 ± 0.08 at = 0.05, and a^/vsed = 0.9 ± 0.1 at </> = 
0.13. All of these values reflect height-averaged values of 
the falling columns, as described above. Ladd et al.'^also 
conducted Lattice Boltzmann simulations of up to 32, 000 
particles sedimenting in a fully bounded box at 0o — 
0.13. They found Oyjvsed — 0.6 for the largest systems. 
All of these sedimentation values are in reasonably good 
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FIG. 16: Inset: Height averaged values ^(2:) vs. particle size a. 

— —1/3 

The dashed line is a fit to ^ = 12a(j)g'". Mainfi gure: velocity 
correlation lengths vs. volume fraction The dashed line 
represents the sedimentation scaling equation, ,^ — llaifiQ , 
for the height averaged correlation length^. 



agreement with our fluidization results shown in Fig. 13 
where near to the middle of the columns, 0.6 < cr^/vp < 
0.9. 



IV. EVALUATION OF THE MUCHA AND 
BRENNER STRATIFICATION MODEL 



The main aim of section IIIII is to describe results for 
the particle dynamics during stable fluidization of a va- 
riety of different sized particle beds. We are not aware of 
other published results or predictions, either experimen- 
tal or theoretical, for the specific forms of the quantities 
we measured, 4>{z), V0(z), (T„(z) or ^(z), in similar flu- 
idized beds. Published results are available for some of 
these properties in similar systems that are not fluidized, 
but rather are sedimenting downwards. Significantly, we 
were able to make a connection between fiuidization and 
sedimentation by showing in Fig. [15] that the local ve- 
locity fiuctuations in a typical sample are very similar in 
either mode. A subsequent comparison confirms this con- 
nection by showing that our average correlation lengths 
(^(z)) and velocity fluctuations in fluidized beds agree 
quite well with published results in sedimenting systems. 
Building upon this link, and in the absence of relevant 
fluidization theories, we therefore flnd it worthwhile to 
turn to theories of particle sedimentation to help explain 
our data. 

The Mucha and Brenner modeP^EH, as outlined in 
the introduction, is a recent theory of sedimentation that 
argues that the value of the local velocity correlation 
lengths ^ are primarily determined by the local degree 
of stratification V(f>, see Eq. ([Tl. Results in support of 
this model, from numerical simulations^, have been con- 
fined to the transition region near the top particle inter- 
face, but whether it can provide a complete description 



of ^ over the entire sedimenting column remains an op^fl 
question. The importance of the MB model, if verified, is 
that it would provide a new insight into the origin of the 
heretofore poorly understood fundamental lengthscale ^ 
that is ubiquitously seen in both fiuidized and sediment- 
ing systems. 

For a complete test of the MB model, in each of 
the four beds we input fits to our data for (/'(z) and 
V(/)(z) into Eq. ([T]) to generate model predictions for 
the velocity correlation lengths l{z)/a. These values are 
then directly compared to our correlation measurements 
^{z)/a. Figure 



14 



shows the results. The solid lines are 
the MB model predictions, and the symbols our data. 
Significantly, there is very good agreement between the- 
ory and experiment for all four fluidized beds. The ad- 
justable fit constant cq, which is not specified by the the- 
ory but is expected to be of order 1, has a single value, 
Co = 0.82, for the 3 scaling samples, a ^ a* . For the non- 
scaling sample, with a ^ a* , the model similarly works 
well, but requires a slightly lower fit value, cq = 0.70. 
The agreement is all the more impressive given that the 
gradients V(/)(z) are highly non-linear, strongly varying, 
functions of height z. The model Eq. ([T]) can also be 
simplified further by noting that the term [0S'(0)]^/^ is 
nearly constant over a wide range of volume fractions, 
[0S'(0)]i/5 ^ 0.52 ± 0.03 for 0.04 < <j} < 0.30. In the 
scaling regime, a 3> a*, therefore 



^(z)/a w O.43[-aV0(z)]-2/5. 



(2) 



The truly remarkable feature about the MB model is 
that it accounts for the characteristic length scale of the 
velocity fluctuations, ^, from which the magnitudes of 
the veloci ty fl uctuations and diffusion coefficients can be 
calculated^, solely from measurements of the stationary 
concentration profiles 4>{z). Its validity corroborates the 
basic assumption of the model, that the velocity fluctu- 
ations arc controlled by the concentration stratiflcation 

V0. 



A. Re-derivation of the MB Model 

Finally, we show that the MB model Eq. ([TJ, ob- 
tained by Mucha and Brenner on the basis of physical ar- 
guments concerning buoyancy matching, can also be de- 
rived from a general advection-diffusion model of steady 
state sedimentation. In this model, as originally pro- 
posed by KynchPS, the local particle concentration varies 
in time as 



d<j)/dt^V[(l)v + DW(j)]. 



(3) 



In our experiments, we achieve stable fluidization so 
that, when averaged over the short fluctuation lifetimes 
T (~ ^/(T„ ~ 5 — 10 s), the concentration profllcs (j){z) 
shown in Fig. [9] are time-independent quantities. Ex- 
pressed mathematically, d{(j))r/dt = 0, where (•••),- rep- 
resents an average over times greater than r. For a stably 



fluidized system, Eq. ^ predicts that the time averaged 
advective, Ja = {4>v)^, and diffusive, Jd — — (DV(f>)^, 
terms are equal at all positions in the particle column. 

To evaluate the advective term, we expand to first 
order, 0(z, t) — (j){z) + S4>{z, t), and v{z, t) — {v{z, t))^ + 
5v{z,t). (For compactness, we will simplify the nota- 
tion below by not writing the explicit z dependence, i.e. 
4i{z) — + (j), but all quantities in general depend upon z). 
Since the mean particle velocity in a fluidized bed van- 
ishes, (f)r — 0, and the time averaged fluctuations are 
zero, {5(j))r = and {5v)r = 0, the advective term re- 
duces to J A = {S4'Sv)^, which is non zero because ve- 
locity and concentration fluctuations are ant i-correlated 
(denser regions fall, lighter regions rise)^'-'. If we ap- 
proximate the fluctuation terms by their rms values, i.e. 
6<j> ^ cr^ = y/{{S(j>)^)^, then 



Ja 



(4) 



To evaluate the diffusive term, J]j, we us e the pub- 
lished result (from sedimentation experiment d^^^^l ) that 
the hydrodynamically induced particle diffusion coeffi- 
cient is related to the fluctuations through D « tr^^, so 
that 



Jd ~ -cr^^Vc?!). 



(5) 



For the stability criteria, J a = Jd, the velocity terms 
ay cancel out, leaving the simple relation Vcj) = — cr^/^, 
that relates the stratification to the concentration fluctu- 
ations. To evaluate tr^, we use the same Poisson model 
as Mucha and Brenner, = y/^i>SJ~4)jaF]~^ . Here, S{4>) 
is the structure factor of random hard spheres needed 
at moderate to high concentrations to account for ex- 
cluded volume effect^. Finally, we solve for the correla- 
tion length ^ and obtain 



(6) 



which is identical to the MB model Eq. ([TJ. 

The new insight provided by this re-derivation of the 
MB model is that it explicitly shows that Eq. (Ill) is con- 
tingent upon steady state behavior, in which d(^r/dt = 
0. This highlights the advantage of testing this model in 
fluidized beds, which produce time invariant concentra- 
tion profiles. In sedimentation on the other hand, several 
studies have reported that the concentr ation pro file con- 
tinues to evolve as the column is falling^ISIlSIiZl^ Conse- 
quently, the profile may or may not reach a steady state 
before the column collapses to the bottom. This deriva- 
tion explicitly shows that in columns with time varying 
concentration profiles, i.e. d(j>{z,t) /dt ^ 0, the advective 
and diffusive particle flux terms would not be balanced, 
and the MB model Eq. ([!]) would not hold true. 



V. DISCUSSION AND CONCLUSIONS 

In this article we present results on the particle dy- 
namics in liquid fluidized beds at low Reynolds number. 



Our primary aim is to examine the behavior of the fluctU 
ating particle velocities and the stationary concentration 
proflles. We chose to do this in moderately concentrated 
samples, all of mean concentration 0o = 0.10. To exam- 
ine how the velocity and concentration profiles depend 
upon system size, we used a single cell geometry, but 
varied the particle size a by a factor of ~ 4.5 so that the 
system size N varies by almost two orders of magnitude. 

For all particle sizes, the measured concentration gra- 
dients Vcj), fluctuation magnitudes ay, and correlation 
lengths ^, all display z— dependent profiles in the particle 
columns. These profiles are found to change in magni- 
tude, but not form, when the particle sizes a are varied 
over a wide range. Specifically, we find a scaling regime 
for particles larger than a critical size, a ^ o*, in which 
the data profiles from different particle sizes collapse to- 
gether when plotted in the forms — aV(/)(z), ^(z)/a, and 
ay{z)/vp. (see Figs. 12 b), 13 and 14 1. These simple 



scalings with particle size demonstrate the universality 
of the particle dynamics and structure in fluidization. 

The new scalings also quantify how the bed proper- 
ties depend upon system size. In our experiments, the 
number of particles fluidized, N = TWH(j)o/{4:/3)TTa^ oc 
1/a^, increases with decreasing particle size. The scaling 
region is found when a 3> a*, or equivalently <C N* , 
where N* « 35 million particles. In this regime, the 
stratiflcation scaling of aV(/)(z) can similarly be writ- 
ten as a N~^^'^Vtf>{z) scaling, explicitly showing that the 
concentration gradient grows with N. Interestingly, in 
the limit of very small systems, the gradients vanish, i.e. 
limjv^oV0(z) = 0, resulting in what would be a uni- 
formly concentrated particle bed. This is an important 
point that shows that it would generally be easier to ob- 
serve concentration gradients in experiments, as opposed 
to simulations, since experimental system sizes are usu- 
ally much larger. On the other hand, Ladd^^ simulated 
the sedimentation of ^ 32000 particles (in a square cell 
geometry) at (po — 0.13 and, while looking across the 
middle of the column, found evidence for an extremely 
small gradient — aVc/) " of order 10~^" . This result agrees 



well with our measurements shown in Fig. 12 where near 
mid-height, -aVcf) w 2 - 3 • 10"^. 

Evidence is also presented that the particle dynam- 
ics in fluidized beds can be closely connected to that seen 
in similar systems that are sedimenting downwards. We 
flrst demonstrate this for a typical sample in Fig 



15 



where the magnitude of the velocity fluctuations show no 
discernible difference between the two modes. To pur- 
sue this further, we apply the same column averaging 
methods to our data that has been used in sedimenta- 
tion experiments to extract a single value for ay and ^ in 
each of our fluidized beds. These values were found to be 
in good agreement with published findings in sedimenta- 
tion, further connecting the two. It is known, however, 
that clear differences do exist between fluidization and 
sedimentation, particulary regarding whether or not the 
systems are in a steady state. Numerous studiej^* ^''*^^ ! 
beginning with GuazzellP^, have shown that the particle 



dynamics (a^, f , Vcfi . . .) in sedimentation can be strongly 
time-dependent as the particles fall, in which case a 
steady state is not achieved. In fluidized systems, on the 
other hand, as the column stays suspended indefinitely, 
a steady state is eventually reached, as we find here. We 
note that in our fiuidization/sedimentation comparisons, 
we compared to sedimentation data thought to be in near 
steady state conditions. 

A further, and central, finding in this paper is the 
verification in Sec. II VI of the Mucha and Brenner model 
for the velocity correlation lengths S.{z). The model pre- 
dicts that ^{z) can be calculated solely from the concen- 
tration profile (j){z) (and its derivative V(j){z)). We mea- 
sured (j){z), and found that the MB model predictions for 
£_{z) from Eq. ^ were in very good agreement with our 
correlation length measurements in all four beds, and at 
all heights within each bed. The physical argument used 
by Mucha and Brenner is that, owing to buoyancy forces, 
concentration gradients limit the maximal size over which 
particle velocities can be correlated, and this leads to the 
relation between ^ and Vcj) given in Eq. ^. In sec. 
IV A we also show an alternative derivation of Eq. ([l]) , 



starting from a general advection-diffusion equation. The 
physical picture we used is seemingly different from the 



MB line of reasoning, yet both yield the same relatiAS 
between ^ and V(j). In our derivation, when the system 
is in a steady state, i.e. d{(j){z,t))r/dt = 0, the diffusive 
flux generated by concentration gradients, D\7(j), is coun- 
terbalanced by an advective flux due to the fluctuations, 
(j^(Ty, which leads to Eq. 

Several key questions still remain to be answered, 
starting with the observation that the MB relation Eq. 
([T| requires a knowledge of the stratification in Vcj) for a 
determination of the correlation length ^, but does not 
allow for an a priori calculation of either quantity. Ad- 
ditionally, while we found there to be an overall system 
size dependence for several properties in steady state sed- 
imentation, we did not vary the different dimensions of 
the cell (TWH) independently, so it is not known whether 
some dimensions of the cell are more influential in this 
regard than others. Future experiments using different 
aspect ratio cells are being planned to help resolve these 
issues. 

We thank Shang Tee and Tony Ladd for stress- 
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ner for many discussions. (corresponding author: 
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